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Abstract 

In hadrons and nuclei at very small x, parton distributions saturate at a scale Qs{x). 
Since the occupation number is large, and Qs{x) >> A^qcd, classical weak coupling 
methods may be used to study this novel regime of non-linear classical fields in 
QCD. In these lectures, we apply these methods to compute structure functions in 
deeply inelastic scattering (DIS) and the energy density of gluons produced in high 
energy nuclear collisions. 
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1 Introduction 



One of more interesting problems in perturbative QCD is the behaviour of parton 
distributions at smaU values of Bjorken x. In deeply inelastic scattering (DIS) 
for instance, for a fixed » Aqq^, the operator product expansion (OPE) 
eventually breaks down at sufficiently small x |jl|. Therefore at asymptotic energies, 
the conventional approaches towards computing observables based on the linear 
DGLAP equations are no longer applicable. Even at current collider energies such 
as those of HERA, where the conventional wisdom is that the DGLAP equations 
successfully describe the data, there is reason to believe that effects due to large 
parton densities are not small. We may be at the threshold of a region where 
non-linear corrections to the evolution equations are large [^, ^. 

In recent years, a non-OPE based effective field theory approach to small x 
physics has been developed by Lipatov and collaborators [^. Their initial efforts 
resulted in an equation known popularly as the BFKL equation [^, which sums 
the leading logarithms of a5log(l/x) in QCD. In marked contrast to the leading 
twist Altarelli-Parisi equations for instance, it sums all twist operators that contain 
the leading logarithms in x. The solutions to the BFKL equation predict a rapidly 
rising gluon density. Such a rapid rise in the gluon density is seen at HERA ||7| but 
it can also arguably be accounted for by the next to leading order (NLO) DGLAP 
equations with appropriate choices of the initial parton densities Q . 

Moreover, the next to leading logarithmic corrections to the BFKL equation 
computed in the above mentioned effective field theory (EFT) approach are very 
large Q. Recently, as Gavin Salam has discussed in his lectures at this school |1C], 
there has been considerable progress in understanding the source of these large 
corrections. The coUinear enhancements from higher orders to the NLO BFKL 
kernel can be resummed, and this results in more stable estimates for the gluon 
anomalous dimensions, and for the hard pomeron. However, there are effects, not 
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included in such analyses, due to multiple pomeron exchange (non-linear QCD 
effects) that may become important at rapidity scales of interest. For instance, 
running coupling effects in the NLO BFKL equation become important at y ~ 

5 /3 

. However, double (hard) pomeron effects will become important for y 
1/as log(l/a5), a scale that is parametrically larger. How to systematically include 
such effects, which are enhanced by large parton densities, is an open question, and 
novel approaches need to be explored. 

An alternative EFT approach to QCD at small x was put forward in a series 
of papers m, H, 0, m. Our approach, is a Wilson renormalization group 
approach (RG) where the fields are those of the fundamental theory but the form 
of the action, at small x, is obtained by integrating out modes at higher values of 



X. This results in a set of non-linear renormalization group equations [15|. If the 



parton densities are not too high, the RG equations can be linearized, and have 



been shown to agree identically with the BFKL and small x DGLAP equations |14]. 
There is much effort underway to explore and make quantitative predictions for the 
non-linear regime beyond |16| ]. 

In these lectures, we will apply the above EFT to discuss two problems: 



Deeply inelastic scattering at small values of Bjorken x |18|, 



high energy hadronic collisions |19, 20, 21]. 



Both are problems which appear extremely difficult to address in an OPE based 
analysis. They simplify in the regime where x is small and momentum transfers are 
large because, a) the gluon field at small x can be treated classically, and b) weak 
coupling methods apply. Why is this so? The reasons for these to apply have been 
discussed at length by Larry McLerran in his lectures |2^ so we will be brief here. 
At small X, one can define a scale //^ which measures the density of gluons per 
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unit transverse area. One has 

^ IdN 

= --r ^ 

a dy 

where a is the hadronic or nuclear cross section of interest. Here y = yo — ln{l / x) , 
and yo is an arbitrarily chosen constant. When this scale satisfies the condition 
/i 3> as/i » ^QCD, the occupation number of gluons in the hadron is large-thereby 
justifying the use of classical methods. Also, the intrinsic momentum pt ~ asfJ- is 
large. Thus the gluon dynamics, while nonperturbative, is both semiclassical and 
weakly coupled. 

Let us now discuss the classical field approach to small x DIS. The gluon field, 
being bosonic, has to be treated non-perturbatively. This is analogous to the strong 
field limit used in Coulomb problems. Fermions, on the other hand, do not develop 
a large expectation value and may be treated perturbatively. In DIS, to lowest 
order m as, the gluon distribution function is determined by knowing the fermionic 
propagator in the classical gluon background field. In general, this propagator must 
be determined to all orders in the classical gluon field as the field is strong. This 
can be done due to the simple structure of the background field. 

We will derive analytic expressions for the current -current correlator in deeply 
inelastic scattering by summing a particular class of all twist operators. These, 
we argue, give the dominant contribution at small x. At high Q^, they reduce 
to the well known expressions for the leading twist structure functions [^]. For 
light quarks at high Q^, it can be shown explicitly that the classical field analysis 
reproduces the DGLAP evolution equations for the quark distributions at small 
[18|. The power of the technique we use to analyze the problem of DIS at small 
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X is that, unlike the OPE, it does not rely on a twist expansion. 

A similar point can be made about the classical approach to high energy 
hadronic scattering. At very high energies, the dominant contribution to parti- 
cle production is from the interaction of the classical "Weizsacker- Williams" (WW) 



gluon fields of tlie two liadrons or nuclei. To lowest order, the picture is that of 



QCD Bremsstrahlung [43|. Soft gluons can be emitted from the valence quark/hard 
gluon lines, or from the 2^1 diagram of two virtual gluons fusing to produce a 
hard gluon. At small x, in the Fock~Schwinger gauge (x~^A~ + x~j4+ = 0) the lat- 
ter WW contribution is the dominant one. The Weizsacker-Williams contribution 



agrees with the QCD Bremsstrahlung result at small x ^ 46 1. 

It is essential to consider the full non-perturbative approach for the following 
reasons. Firstly, the classical gluon radiation computed perturbatively is infrared 
singular, and has to be cut-off at some scale kt ~ asA*- This problem also arises in 
mini-jet calculations where at high energies results are shown to be rather sensitive 
to the cut-off . Secondly, the non-perturbative approach is crucial to a study of 
the space-time evolution of the system. In particular, the possible thermalization 
of the system, as well as the relevant time scales for thermalization, are strongly 
influenced by the non-linearities that arise in the non-perturbative approach [p^ j. 

In these lectures, we will discuss results from real time simulations of the full, 
non-perturbative, evolution of classical non-Abelian WW fields. The fields are gen- 
erated by sources of color charge (representing the valence partons in each of the 
hadrons or nuclei) moving along the two light cones. For each p configuration, one 
solves Hamilton's equation numerically to obtain the real time behavior of the gauge 
fields in the forward light cone. The Hamiltonian is the Kogut-Susskind Hamilto- 
nian in 2-|-l-dimensions coupled to an adjoint scalar field. The initial conditions 
for the evolution are provided by the non-Abelian Weizsacker-Williams fields for 
the nuclei before the collisions. 

To compute observables, one has to average over all the p configurations in 
each of the two nuclei. In general, these are averaged with a statistical weight 
exp [— where F[p] is a functional over the color charge density p of the higher 
X modes. The functional F[p] obeys the non-linear renormalization group equation 
that was mentioned in the preceding discussion. If one considers collisions of large 
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nuclei, the weight simplifies to a Gaussian one, and one can replace 

where yf^ag the rapidities corresponding to the fragmentation regions of the two 
nuclei. 

Our approach is limited because it is classical. However, if the effective action 
approach captures the essential physics of the small x modes of interest, then in the 
spirit of the Wilson renormalization group, quantum information from the large x 
modes (above the rapidity of interest) is contained in the parameter fj,'^ discussed 
above, which grows rapidly as one goes to smaller x's. In principle, this information 
can be included in the classical lattice simulations. 

The plan of these lectures is as follows. We will begin in section 2 by reviewing 
the effective action for the small x modes in QCD. We also discuss the classical 
saddle point solutions of this effective action. In section 3, we will discuss how one 
computes quark production in the classical gluon background field. At small x, 
this gives the dominant contribution to the structure functions measured in deeply 
inelastic scattering. Structure functions are computed in section 4. At high Q^, 
our results reproduce the small x DGLAP results. For smaller values of Q^, for 
Gaussian sources, one obtains the Glauber result for the structure functions in 
agreement with previous derivations in the nuclear rest frame. Subsequent sections 
concern gluon production in high energy scattering of (in particular) large nuclei. 
The classical approach to the two-nucleus problem is discussed in section 5. It is 
very hard to solve for the non-perturbative dynamics analytically. It has not yet 
been done. Instead, we derive a numerical algorithm which captures the essential 
physics of the two-nucleus problem. Results from our numerical simulations are 
discussed in section 6. Section 7 summarizes the material contained in the lectures 
and outlines directions of future research. 
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2 Effective Field Theory for Small x Partons in QCD 



In the infinite momentum frame (IMF) P+ oo, the effective action for the soft 
modes of the gluon field with longitudinal momenta A;^ << (or equivalently 
X = k^/P^ << 1) can be written in light cone gauge A'^ = as 

Seff = - j d^x-^G^G^'^^^ + j^ j d^xtdx-p^{xux-)l^{T^W^^^^[A-]{x-,xt)) 
+ i j d^xtdx-Flp^ixux-)]. (3) 

Above, is the gluon field strength tensor, t"" are the SU{Nc) matrices in the 
adjoint representation and W is the path ordered exponential in the x^ direction 
in the adjoint representation of SU{Nc), 



W^oo,oo[-A ]{x ,xt) = Pexp 



The action is a gauge invariant form |14] of the action that was proposed in Ref. |11]. 
One can write an alternative gauge invariant form of the action but the results are 
the same for the problem of interest. 

The effective action considered here is valid in a limited range of P^ « A^, 
where A^ is an ultraviolet cutoff in the plus component of the momentum. The 
degrees of freedom at higher values of P^ have been integrated out. Their effect 
is to generate the second and third terms in the action. The first term is the 
usual field strength piece of the QCD action and describes the dynamics of the 
wee partons. The second term is the coupling of the wee partons to the hard color 
charges at higher rapidities, with x values corresponding to values of P'^ > A"*". 
When expanded to first order in yl~, this term gives the familiar J ■ A coupling for 
Abelian classical fields. The last term in the effective action is imaginary. It can be 
thought of as a statistical weight resulting from integrating out the higher rapidity 
modes in the original QCD action. Expectation values of gluonic operators 0{A) 
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are then defined as 

. _ /[^p]exp(-F[p])/[dA]0(A)exp(iS[p,^]) 

/[dp]exp(-FM)/[dA]exp(.5[p,^]) ' 

where A] corresponds to the first two terms in Eq. ^. 

In the IMF, only the J+ component of the current is large (the other components 
being suppressed by 1/P+). The longer wavelength wee partons do not resolve the 
higher rapidity parton sources to within 1/P+ and, for all practical purposes, one 
may write 

p\xt,x-) ^ p\xt)b{x-) . (6) 



In Ref. 1 11] a Gaussian form for the action 

d'x.^p^p^ (7) 

was proposed, where p} was the average color charge squared per unit area of the 
sources at higher rapidities. For large nuclei A >> 1 it was shown that 

1 iVo 



,,2 ^ ^'g 
p = ~ 



A^/ye fm~^. (8) 



This result was independently confirmed in a model constructed in the nuclear rest 



frame |23]. If we include the contribution of gluons which have been integrated out 



by the renormalization group technique, one finds that |^T[ 

p^ = ^(^ + ^^] (9) 

Here Nq is the total number of quarks with x above the cutoff; Nq = J2i fx dx'qi{x') 
where the sum is over different flavors, spins, quarks and antiquarks. For gluons, 
we also have Ng = dx'g{x'). The value of ttR'^ is well defined for a large nucleus. 
For a smaller hadron, we must take it to be a, the total cross section for hadronic 
interactions at an energy corresponding to the cutoff. This quantity will become 
better defined for a hadron in the renormalization group analysis. 



The above equation for /x^ is subtle because, implicitly, on the right hand side, 
there is a dependence on /x through the structure functions themselves. This is the 
scale at which they must be evaluated. Calculating /i therefore involves solving an 
implicit equation. Note that because the gluon distribution function rises rapidly 
at small x, the value of /i grows as x decreases. At some critical x, the indications 
are that the parton distributions saturate. Thus there may be a critical line in the 
x-Q'^ plane corresponding to parton saturation. This is an important point and we 
will return to it later. 

The Gaussian form of the functional F[p] is reasonable when the color charges 
at higher rapidity are uncorrelated and are random sources of color charge. This 
is true for instance in a very large nucleus. It is also true if we study the Fock 
space distribution functions or deep inelastic structure functions at a transverse 
momentum scale which is larger than an intrinsic scale set by a^/x. In this equation 
as is evaluated at the scale /x. At smaller transverse momenta scales, one must do 
a complete renormalization group analysis to determine F[p]. For heavy quarks, in 
DIS, the Gaussian analysis should be adequate. 

In Ref . , it was shown that a Wilson renormalization group procedure could 
be applied to derive a non-linear renormalization group equation for F[p]. In the 
limit of weak fields, the renormalization group equation can be linearized, and can 
be shown to be none other than the BFKL equation discussed previously. The fact 
that this limit can be obtained in a simple and elegant way suggests the power 
of this approach, and the importance of further studying the non-linear region of 
strong classical fields. We will not discuss the RG procedure here but will refer the 



reader to the relevant papers, and to Larry McLerran's lectures [22 

The effective action in Eq. ^ has a remarkable saddle point solution |11, 13, ^3f| . 
It is equivalent to solving the Yang-Mills equations 

D^G^'^ = r6''+ , (10) 
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in the presence of the source J"*"'"^ = p^{xt, x~). Here we will allow the source to 
be smeared out in x~ as this is useful in the renormalization group analysis. It is 
also useful for intuitively understanding the nature of the field. One finds a solution 
where = and 

A' = — Vd'V^ , (11) 

(i = 1, 2) is a pure gauge field in the two transverse dimensions which satisfies the 
equation 

Di—=gpiy,x±). (12) 

Here Di is the covariant derivative di + VdiV^ and y = yo + log(x~/xg ) is the 
space-time rapidity and i/q is the space-time rapidity of the hard partons in the 
fragmentation region. At small x, we will use the space-time and momentum space 



notions of rapidity interchangeably [17|. The momentum space rapidity is defined 
to be y = 2/0 — ln{l/x). The solution of the above equation is 



1 



The classical nuclear gluon distribution function is computed by averaging over 
the product of the classical fields in Eq. |l^ at two space-time points with the weight 
F[p] [^]. For a Gaussian source, one obtains 

where Ci;' is the Casimir in the fundamental representation and a is the nuclear 
cross-section|^. For large xt (but smaller that 1/Aqcd, the distribution falls like a 

power law l/x|-and has a 1/as dependence! For very small xt, the behavior is the 
^ Above, we have re- written the expression for the gluon distribution in Ref. using the 

leading log gluon distribution to replace fi^ and log(a;tAQCD) with the gluon distribution xG{x, ■\) 
at the scale l/x^. 
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perturbative distribution log{xtAQCD)- The scale which determines the cross-over 



from a logarithmic to a power law distribution is, following Mueller's notation |24], 
the saturation scale Qs- Setting xt = 1/Qs and the argument of the exponential 
above to unity, one obtains the relation, 

Ql = ^-^xG{x,Ql), (15) 

which, for a particular x, can be solved self~consistently to determine Qs- 

Because of the sharp cut-off in co-ordinate space, the momentum space dis- 
tribution is not well defined. A smooth Fourier transform has been defined, on 
physical grounds, by Lam and Mahlon by requiring that the charge in light cone 
gauge / d?Xtp{xt) vanish [^] for each p configuration. 

3 Quark production in the classical gluon background 
field 

In this section, we will compute the correlator of electromagnetic currents in the 
classical gluon background field. In deeply inelastic electroproduction, the hadron 
tensor can be expressed in terms of the forward Compton scattering amplitude r^,y 



by the relation [26| 



W^'^iq^, P-q) = 2Disc T^^iq^.P ■ q) = ^Im / (fx exp{iq ■ x) 

X < P|r(J^(x)J^(0))|P >, (16) 

where "T" denotes time ordered product, = 'ipj^ip is the hadron electromagnetic 
current and "Disc" denotes the discontinuity of T^u along its branch cuts in the 
variable P • q. Also, — > oo is the momentum transfer squared of the virtual 
photon^ and P is the momentum of the target. In the IMF, P~^ — > oo is the 
^ Note that in our metric convention, a space-like photon has = > 0. 
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only large component of the momentum. The fermion state above is normalized as 
< P \ P' >= {27r fE/m6^^\P - P') where m is the mass of the target hadron. 
This definition of W'^'^ and normalization of the state is traditional. In the end, 
all factors of m cancel from the definition of quantities of physical interest. (The 
normalization we will use in this paper for quark and lepton states will have E/m 
replaced by 2P+.) 

We now generalize our definition of W'^'^ to a source which has a position 
dependence. We obtain 

W*'^(g^ P • q) = ^a—lm [ d^xdX'e''^'' < T {J^iX' + x/2)r{X- - x/2)) > . 

^TT 71% J 

(17) 

To see this, first note that we can define < O >=< P\0\P>/<P\P> where 
O is any operator. Prom the discussion above, the expectation value < P \ P >= 
(27r)'^£'/m 5(^^(0) = {27r f E/m V. Here we shall take the spatial volume V to be 
(7 times an integral over the longitudinal extent of the state. The variable X^ is a 
center of mass coordinate and x~ is the relative longitudinal position. The above 
definition of W^'^ is Lorentz covariant. The integration over X~ is required since 
we must include all of the contributions from quarks at all X~ to the distribution 
function. In the external source language, the variable P'^ can be taken to be the 
longitudinal momentum corresponding to the fragmentation region. 

The expectation value is straightforward to compute in the limit where the 
gluon field is treated as a classical background field. If we write 



< 



TiJ'^{x)r{y)) >=< T ['iljixh^^lj{xmy)Yip{y)) > , (18) 
then when the background field is classical, we obtain 

< r(J'^(x)J'^(y)) >= TTirSA{x))TT{YSAiy)) + Tiir SAix,y)YSA{y,x)) . (19) 
In this expression, SA{x,y) is the Green's function for the fermion field in the 
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external field A 



SA{x,y) = -i<i^{x)ij{y)>A . (20) 



The first term on the right hand side of Eq. |19| is a tadpole contribution without 
an imaginary part. It therefore does not contribute to W^^ . We find then that 

W^''iq^,p-q) = — a— Im [ dX-d^xe'i'' (Tr(j'' Sa(X~ +x/2,X- -x/2) 
271 m J ^ 

X ^''Sa{X- -x/2,X- + x/2))). (21) 



The expression we derived above for W^'^ is entirely general and makes no 
reference to the operator product expansion. In particular, it is relevant at the small 
X values and moderate where the operator product expansion is not reliable [|^. 
At sufficiently high though (and for massless quarks) it should agree with the 
usual leading twist computation of the structure functions. 

We can derive an expression for the sea quark Fock distribution in terms of the 
propagator in light cone quantization |2?]. One obtains 



dN 2i 



d^k {2TTy 



J d^x d^y e-'^-^'^-y'^ Tr [7+5a(x, y)] (22) 



where the fermion propagator Sa{x, y) is defined as in Eq. 20. In a nice pedagogical 



paper, (see Ref. [28| and references within), Jaffe has shown that the Fock space 
distribution function can be simply related to the leading twist structure function 
F2 by the relation 

F2(x,Q2)= dk^x^^. (23) 

Actually, Jaffe's expression is defined as the sum of the quark and anti-quark dis- 
tributions. At small X, these are identical and the resulting factor of 2 is already 
included in our definition of the light cone quark distribution function. At high q^, 
our general (all twist) result for F2 agrees with the leading twist result derived using 
Eq. II 0. 
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Clearly, to compute W^'^, we first need to compute the fermion Green's function 
in the classical background field. The field strength carried by these classical gluons 
is highly singular, being peaked about the source (corresponding to the parton 
current at x values larger than those in the field) localized at x~ = 0. Away from 
the source, the field strengths are zero and the gluon fields are pure gauges on both 
sides of x~ = (see Eq. ^). The fermion wavefunction is obtained by solving the 
Dirac equation in the background field on either side of the source and matching the 
solutions across the discontinuity at x~ = 0. Once the eigenfunctions are known, 
the fermion propagator can be constructed, in the standard fashion, by writing 



after identifying = [q^ + — X)/2q^. We will not discuss the details of the 



derivation here but refer the interested reader to Ref. 
Define 

G{xt,x-) = e{-x-) + e{x-)V{xt) , (25) 

a gauge transformation matrix that transforms the gluon field at hand to a singular 
field with the only non-zero component, A'^ = 5^^a{xt)- Our result then is that 
the fermion propagator in the background field has the form |29, 18 1 



SA{x,y) = G{x)So{x-y)G\y)-i j d^zG{x)\^e{x~)e{-y-){vHzt)-l)- 
e{-x-)e{y'){V{zt) - l)^G\y)So{x - z)-i-5{z-)So{z - y) . (26) 
with the free fermion Green's function 

SJx -v)- f c-g-fa-^) (^-^) (27) 
•^"^"^ y'~J (27r)4" q2+M^-ie ■ 

Recall that V{xt) is the gauge transformation matrix in the fundamental represen- 
tation and that the classical solution = V{xt)d^V^ {xt)/{—ig). This very simple 
form of the propagator is useful in the manipulations below. 
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In fact, since the current-current correlation function is explicitly gauge invari- 
ant, we may use the singular gauge form of the propagator |^, for computing 
the current-current correlation function 

5r'(x,y) = So{x-y)-iJd'zi^e{x-)9{-y-){VHzt)-l)- 

e{-x')e{y-){V{zt) - l)|so(x - z)tKz-)S^{z - y) . (28) 

A diagrammatic representation of the form of the propagator above is shown in 
Fig. 1 In the expressions below for W^^ we will drop the superscript sing and 
simply use the singular gauge expression for the propagator. 




Figure 1: Diagrammatic representation of the propagator in Eq. 48. 



Our result for the fermion propagator in the classical background field was 
obtained for a J-function source in the x~ direction. This assumption was motivated 
by the observation that small x modes with wavelengths greater than 1 /P^ perceive 
a source which is a (5-function in x~ . The propagator above can also be derived for 
the general case where the source has a dependence on x~ . The gauge transforms 
above are transformed from V{xt) V{xt, x~), to path ordered exponentials, where 
V{xt, x^) is given by Eq. [l^. Our result for the propagator is obtained as a smooth 
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limit of Ax^ = » x~(= 1/P^). Therefore our form for the propagator is 

the correct one provided we interpret the 0-functions and 5-functions in x~ to be 
so only for distances of interest greater than , the scale of the classical source. 

We are now in a position to calculate the current-current correlator. This 
calculation is accurate to lowest order in as but to all orders in asA*- Before we go 
ahead with the computation, we will discuss briefly the averaging procedure over 
the labels of color charges at rapidities higher than those of interest. This is required 
if we are to compute gauge invariant observables. 

If we average the Green's function in Eq. over all possible values of the color 
labels corresponding to the partons at higher rapidities, we can employ the following 
definitions for future reference. Defining 

< Tr {y{xt)V\yt)) >p= ^{xt - yt) , (29) 

we see that 7(0) = 1, which follows from the unitarity of the matrices V. Now 
defining the Fourier transform^ 

^{pt) = J (fxt e-*P*^* Uxt) - 1] , (30) 

we have the sum rule 



I 



7(pt)=0. (31) 



(2vr)2 

The function ^{pt) will appear frequently in our future discussions and as we shall 
see, can be related to the gluon density at small x. 

We will now use the fermion Green's function in Eq. ^ to derive an explicit 
result for the hadronic tensor W'^'^ . In the following section, we will compute the 
structure functions Fi and F2. As previously, we define 

W^'^iq, P, X-) =lm f (fz e'"-' < r( J''(X- + ^)r{X- - J) > , (32) 



^We define the Fourier transform in this way because it corresponds to only the connected pieces 
in the correlator. 
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where "Imaginary" stands for the discontinuity in q . Then 

W'"'(q,P) = —a— I dX-W^(q,P,X-) = —aP+lm f dX' f d^ze''^-' 
2tt m J 2tt J J 

X Tr(sAjx~ + ^,X--^)YSA^,(x--^,X~ + ^)r) ■ (33) 



2 2 J V 2 2, 

The only terms in the propagator that contribute to the above are the 0(ibx~)0(=Fy ) 
pieces. Using our representation for the propagator in Eq. after considerable 
manipulations, we can write W^'^ as 

M//^-r ^:^T f "^'P -fi.^ 



X Tr ■ 



(M - ]/)7-(M - f)Y{M - l'/)T{M - p'/H 



\ (p2 + M2 - ie) (/2 + M2 - ie) (/'^ + M2 - ie) (p'^ + m2 - 



(34) 



where l=p — k, l' = l — q, p'=p — q and /c~ = 0. Correspondingly, we can write 
W^'^ as the imaginary part of the diagram shown in Fig. 2. 

For the DIS case, q^ > (see footnote 1), we can cut the above diagram only in 
the two ways shown in Fig. 3 (the diagram where both insertions from the external 
field are on the same side of the cut is forbidden by the kinematics) . 

Also interestingly, the contribution to W^'^ can be represented solely by the 
diagram in Fig. 4 and not, as is usually the case, from the sum of this diagram and 
the standard box diagram. This is because in our representation of the propagator 
multiple insertions from the external field on a quark line can be summarized into 
a single insertion. Eq. |2^ makes this point clear. 

Applying the Landau-Cutkosky rule, shifting p p + k, and changing variables 



appropriately, Eq. 34 can be written as 



X (27rf6iip + kf+M^)6{{p-q)^ + M^)) ^ ^ 



pi + M'^ {p+k- qY+M"^ ' 

(35) 
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Figure 2: Polarization tensor with arbitrary number of insertions from the classical 
background field. Wavy lines are photon lines, the solid circle denotes the fermion 
loop and the dashed lines are the insertions from the background field (see Fig. 1). 
The imaginary part of this diagram gives W^'^. 



where above the trace is represented by 

^f^u = Trj (M - - f^)-f- (Af - i/h^iM - V + g')7" (M - j/ - ^ + (^)Y + ^ v 



In Ref. [18 1, the integral over p in Eq. 34 was performed before the p'^ integral. 
Here, we instead perform the p"*" integral first. Further, defining z = p^ /q~ , we 



note that the ^-function and 5-function constraints in Eq. 34 restrict < z < 1. 



^Kinematic note: the observant reader will notice we have put — Q here. Since we are 
working in the infinite momentum frame, the hadron has only one large momentum component, 
P+. The rest are put to zero. For the photon, we choose a left moving frame such that q = \q^\ 
and g+ = 0. Then, q^ = qf > 0, P ■ q = -P+q- and Xbj = ~q^/{2P ■ q) = gf/(2P+g"). Since in 
the infinite momentum frame < xsj < 1, this gives q~ > 0. 
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p-q 



p — q + k 




Figure 3: Cut diagrams corresponding to the imaginary part of W^'^. 



This simphfies the result in Ref. []18[ considerably. Performing the p~^ integral, one 
obtains 



1 



aP+q-Nr 1 



1 



dz 



(fipt dPkt .M^"' 



(27r)2 (27r)^ 
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M2 + z(l-z)g2 Ml^^_^ + z{l-z)q^ 



2 ' 



(37) 



where = p^ + . Similarly Mp_^^,_g = {pt + h - qtf + M^. The above is the 
final result of this section, and will be used below in computing structure functions. 

4 Structure Functions at Small x 

The hadronic tensor W^^ can be decomposed in terms of the structure functions 

q^'{P■ 



Fi and F2 as |26 
mW" = -{ c 



q^P ■ q\ F2 



(38) 



where P^ is the four-momentum of the hadron or nucleus and P^ = m? ~ 
(<< q^). In the infinite momentum frame, we have P~^ — > 00 and P~ , Pt ~ 0. 
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Figure 4: In the singular gauge representation for the propagator (see Eq. 28 and 
Fig. 1), the current-current correlator (imaginary part of LHS) is equivalent to the 
imaginary part of RHS. 

The above equation can be inverted to obtain expressions for Fi and F2 in terms of 
components of W^'^ . Since in our kinematics = (see footnote 3 for a kinematic 
note) we have 

It is useful to verify explicitly that our expression for W^'^ derived in an external 



field can be written in the form of Eqn. |38|. Recall that W'^'^ can be written in 
Lorentz covariant form by using the vector = 5^~^. Using n-j = —7" in Eqn. 
we see that W^'^ is a Lorentz covariant function of the only vectors in the problem- 



q^^ and n^. Identifying = P^/P^ in Eqn. 38, we see that these forms are in 



complete agreement. Note that all factors of m disappear from Fi and F2 by the 



explict forms of Eqns. 38 and 34 



We also see that the structure functions can only be functions of and n • q 
by Lorentz invariance. We can therefore take q^ = for the purpose of computing 

20 
1 



Fi and F2. 

To compute TF"*"^ and W , we need to know the the traces Af++ and M , 
respectively in Eq. ^ We can compute them exphcitly and the results can be 
represented compactly as 

^ M++ = I (m2m2^,_^ + M2+,M2_^ - q^k^) , (40) 
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and 



M— =32{p-f{p- - q-f . (41) 

From the relations above of Fi and F2 to W~^~^ and W , we obtain from Eq. ^ 
the following general results for the structure functions for arbitrary values of Q'^, 
and the intrinsic scale 

-^2= = ^^^^ dz dxtXt {l-'^{xt,y)) 

X 'Kl{xtA) (4^2(1 _ zfQ^ + M2) + Kl{xtA)A^ (z^ + (1 - zf)\ (42) 

Here = Q'^z{\ — z) + and i^o,i are the modified Bessel functions. For 
simplicity, we have ignored the impact parameter dependence of 7-and replaced 
the integral over impact parameter by the transverse area a. For the same reason, 
we ignore the sum over the charge squared of the quark flavors. Both of these must 
of course be included in numerical computations. The first (second) term in the 
bracket above is the probability for a longitudinally (transversely) polarized photon 
to split into a qq pair. Ignoring target mass corrections which are negligible at small 

Fl = F2-2xFi 

2^AT rl 



dz / ^'^^^^ dxt xt (1 - 7(xt, y)) z\l - zfQ^ Kl{xtA) . (43) 



'which is implicitly contained in the function 7(fet) in Eq. 
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For a Gaussian source (see Ref. [18] and footnote 1), 

l{xt, y) = exp ( ^) ] , (44) 



where the scale is set by the transverse separation xt between the quark and the 
anti-quark. 

The equation for F2 with the Gaussian source is the weU known Glauber expres- 



sion [32 1 usually derived in the rest frame of the nucleus. It is heartening to see that 
the formahsm of Ref. for structure functions in the infinite momentum frame 
reproduces it. For large , it reduces to the standard small x DGLAP expression^ 
while at small it goes to zero as QHog{Q'^). One then recovers, qualitatively, 
the shape of the famous Caldwell plot for dF2 / d\og{Q'^) measured at HERA p^ . 
Similar forms were used by several authors to understand the recent data [^]. 
One obtains from the above equation for in a manner analogous to Eq. H, 



the quark saturation scale Qf by replacing Cp — > Ca in Eq. 15. The relative size 
of the two saturation scales, glue to quark, is therefore determined simply by the 
ratio of the two Casimirs, Ca/Cf- 

What about quantum corrections to the above quark and gluon distributions? 
At the one loop level, one gets log(l/x) corrections to the Weizsacker-Williams 



distribution |35, 12, 24]. However, Mueller has argued recently that beyond the one 
loop level, the distribution has the same form as the as the above classical gluon 
distribution. What does change due to small x evolution is the x dependence of the 
saturation scale W%. Recently, there have been detailed studies by Kovchegov, and 



by Levin and Tuchin, of parton evolution in the non-linear region p6| , 37]. Their 
results appear to confirm the intuitive picture of Mueller. 

As — > 00, we find remarkably that the integral on the RHS of Eq. 



^In Ref. [Q, it was shown explicitly that our general expression for F2 formally reduces to the 
leading twist expression obtained from Eq. E^. 
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vanishes and it reduces to 

F, = £^ (45) 

The above is the well known Callan-Gross relation. The reader may note above that 
the deviation from the Callan-Gross relation vanishes as a power law as — > oo. 



On the other hand, it is well known in QCD [P9| , 40 1 that the violations of the 
Callan-Gross relation only disappear logarithmically as — > co. The apparent 
contradiction is resolved by one realizing that the logarithmic violations at large 
in QCD come from diagrams where the sea quark emits a gluon (thereby violating 
Feynman's parton model helicity argument). These diagrams are of higher order in 
our picture and are therefore not included. In fact, the deviations from the Callan- 
Gross relation of the sort discussed above (at small x) should die off faster than 
logarithmically at very large because for sufficiently large q^, the violations of 
the Callan-Gross relation should come from precisely the diagrams not included 
here. At moderate q^ however, the contributions we have discussed above should 
be important. 

5 The non— Abelian Weizsacker— Williams approach to 
high energy nuclear collisions 

In nuclear collisions at very high energies, the hard valence parton modes in each 
of the nuclei act as highly Lorentz contracted, static sources of color charge for the 
wee parton, Weizsacker-Williams modes in the nuclei. The sources are described 
by the current 

r^'^in) = 5-+gp'{{n)5{x-) + 5^- gp'i{rt)5{x'^) , (46) 

where pi (p2) correspond to the color charge densities of the hard modes in nucleus 1 
(nucleus 2) respectively. The classical field of two nuclei is described by the solution 
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of the Yang-Mills equations in the presence of the light cone sources: 



(47) 



Gluon distributions are simply related to the Fourier transform A^{kt) of the 
solution to the above equation by < Af{kt)Af{kt) >p. The averaging over the 
classical charge distributions is defined by 



(48) 



and is performed independently for each nucleus with equal Gaussian weight /i^. Of 
course, this is only true for identical nuclei. 

Before the nuclei collide (t < 0), a solution of the equations of motion is 



A^ = ; A' = 9{x-)e{-x+)a\{rt) + e{x+)9{-x-)a2irt) . 



(49) 



where a'^g{rt) {q = 1,2 denote the labels of the nuclei) are pure gauge fields of the 
two nuclei before the collision and have the form described in Eq. The above 
expression suggests that for t < the solution is simply the sum of two disconnected 
pure gauges. 

For t > the solution is no longer pure gauge. Working in the Schwinger gauge 



A'^ = x~^A + X A'^ = the authors of Ref. |42] found that with the ansatz 
A^ = ±x^aiT, n) ; A' = ai_{T, n) , 



(50) 



where r = V2x^x , Eq. ^ could be written in the simpler form 

1 



:drr'dra+[D„ 



D\a 



0. 



1 



[Di,dra^^ + igT[a,dra] = 0. 



1 



drTdra'± — igr [a, 



D\a 



]-[D^,F^'] = 0. 



(51) 



The above equations of motion are independent of r/-the gauge fields in the forward 
light cone are therefore only functions of r and rj and are explicitly boost invariant. 
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This result follows from the assumption that the sources of color charge are delta 
functions on the light cone. Of course this is not true in general. However, we 
are interested in the region of central rapidity, about one unit of rapidity around 
rj = 0. The boost invariance assumption should be Ok in this region. Also note 
that boost invariance is not assumed when solving for the fields of the nuclei before 
the collision. 

The initial conditions for the fields a{T,rt) and at r = are obtained by 
matching the equations of motion (Eq. |47| ) at the point = and along the 
boundaries = 0,x~ > and x~ = 0,x~^ > 0. Remarkably, there exist a set of 
non-singular initial conditions for the smooth evolution of the classical fields in the 
forward light cone. These can be written in terms of the fields of each of the nuclei 
before the collision (t < 0) as follows, 

Q;i|r=o = + "2 ; «It=o = -^["1,02] • (52) 



Gyulassy and McLerran have shown [41| that even when the fields a\ 2 before the 



collision are smeared out in rapidity to properly account for singular contact terms in 
the equations of motion the above boundary conditions remain unchanged. Further, 
the only condition on the derivatives of the fields that would lead to regular solutions 
are dra\r=o, draWT=o = 0. 



In Ref. 1 42], perturbative solutions (for small p) were found to order by 
expanding the initial conditions and the fields in powers p (or equivalently, in powers 
of asp/h) We will not discuss the details of the perturbative solution but wish to 
refer the reader to the original papers. 

Perturbatively, at late times, the fields in the forward light cone can be expanded 
out in plane waves. The energy distribution in a transverse box of size R and 
longitudinal extent dz can be computed by summing over the energy of the modes 
in the box with the occupation number of the modes given by the mode functions 
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ai{kt). We have then 

The multiphcity distribution of classical gluons is defined as dE / dyd? kt / u) . After 
performing the averaging over the Gaussian sources, the number distribution of 
classical gluons is 



dyd^h (27r)4 kt 



.n-u ^4 -Lih, A) , (54) 



where L(kt,X) is an infrared divergent function at the scale A. This result agrees 
with the quantum bremsstrahlung formula of Gunion and Bertsch |43] and with 
several later works @ ||, |6|. 

The function L{kt,X) arises from long range color correlations that are cut- 



off either by a nuclear form factor (as in Refs. |43, by dynamical screening 

effects |5^ or in the classical Yang-Mills case of Ref. non-linearities that 
become large at the scale kt ~ asfJ-- In the classical case, L{kt,X) = log(A;|/A^), 
where A = as/U. The formalism used in all these derivations breaks down at small 
momenta and one cannot distinguish between the different parametrizations of the 
nuclear form factors. However, at sufficiently high energies, the behaviour of L(kt, A) 
in the infrared is given by higher order (in asfJ-fkt) non-linear terms in the classical 
effective theory. We hope to understand in the near future how non-perturbative 
effects in the classical effective theory dynamically change the gluon distributions 
at small transverse momenta. 

While the Yang-Mills equations discussed above can be solved perturbatively 
in the limit asfJ- <C kt, it is unlikely that a simple analytical solution exists for 
Eq. ^ in general. The classical solutions have to be determined numerically for 
t > 0. The straightforward procedure would be to discretize Eq. |4^ but it will be 
more convenient for our purposes to construct the lattice Hamiltonian and obtain 
the lattice equations of motion from Hamilton's equations. 
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Let us first consider the continuum Hamiltonian |1S]. In the appropriate (r, rj, xt) 
co-ordinates, the metric is diagonal with g'^'^ = —g^^ = —g^^ = 1 and g^^ = — 
After a little algebra, the Hamiltonian can be written as [^] 

H = tJ dr]d\ jipV + + ^Fr,rFr,r + ^F^yF.y + f + f^} • (55) 

Here we have adopted the gauge A"^ = 0. Also, = ^drA^j and p'" = rdrAr are 
the conjugate momenta. 

Consider the field strength F^/r in the above Hamiltonian. If we assume approx- 
imate boost invariance, or 



Ar{T,r],ft) ~ Ar{T,ft); Ar,{T,r],i%) ~ ^{T,ft), 



we obtain 



(56) 



(57) 



where = dr — igAr is the covariant derivative. Further, if we express j^'^ in 
terms of the defined in Eq. 46 we obtain the enormously simplifying result that 



3 



r],r 



for r > 0. Due to boost invariance, our effective Hamiltonian acts in 



2+1-dimensions. It is possible to relax this assumption, but then the numerical 
simulations are more complicated. 

We now consider the equivalent lattice action and Hamiltonian. The appropri- 
ate action is derived starting from the Minkowski Wilson action in the discretized 
4-space and taking the naive continuum limit in the longitudinal and time direc- 
tions. Replacing with / dzdt in the Minkowski Wilson action, we then have 
for the 2+1-dimensional action 



S = 
where 



J dzdtY, 



_i_TrFi=, + i-!RTr(Ma-M.J 



1 ^TtU± 

Nr. 



(58) 



3,n 



(59) 
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and similarly for M^jn- 

The equation of motion for a field is obtained by varying S with respect to 
that field. For the longitudinal fields At^z the variation has the usual meaning 
of a partial derivative. For transverse link matrices U± the variation amounts to a 
covariant derivative. Just as in the continuum case, the lattice initial conditions can 



be determined from the lattice action in Eq. 58. One obtains the lattice equations 
of motion in the four light cone regions and then determines non-singular initial 
conditions by matching at r = the coefficients of the most singular terms in the 
equations of motion. 

On the lattice, the initial conditions are the constraints on the longitudinal 
gauge potential and the transverse link matrices U± at r = 0. The longitudinal 



gauge potentials can be written as in the continuum case (see Eq. 50) as 



A± = ±x^e{x+)e{x-)a{T, xt) . (60) 

The transverse link matrices are, for each nucleus, pure gauges before the collision. 
This fact is reflected by writing 

= e{-x+)9{-x-)i + 9{x+)eix-)u{T) + 9{-x+)e{x-)Ui + e{x+)e{-x-)U2 , (6i) 

where Ui^2 are pure gauge. The pure gauges are defined on the lattice as follows. 
To each lattice site j we assign two SU(A'^c) matrices Vij and V2j- Each of these 
two defines a pure gauge lattice gauge configuration with the link variables C/J^ = 
VqjVqj+n whcrc q = 1,2 labels the two nuclei. As in the continuum, the gauge 
transformation matrices Vqj are determined by the color charge distribution pqj of 
the nuclei, normally distributed with the standard deviation /i^: 



(62) 

Parametrizing Vqj as exp(iAp with Hermitean traceless Aj, we then obtain Aj by 
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solving the lattice Poisson equation 

E ( Al+n + Al-n - 2 A]) = P,, ■ (63) 



It is easy to verify that the correct continuum solution (Eqs. ^ and 50) for the 
transverse fields is recovered by taking the formal continuum limit of Eq. 

The equation of motion for f7_L, contains, upon substitution of U±_ from (^Tj) and 
from dH), singular terms containing the product 5{x )5{x'^). These originate 
in the double-derivative contributions 5RTrC/j^5-|-9_C/j_ in the action, when both 
derivative operators act on the step functions. Since the coefficient in front of 
6{x^)6{x^) must vanish in order to satisfy the equation of motion, a matching 
relation between U± and Ui^2 is obtained. 

Trcj^, [{Ui + U2){l + Ui) - h.c] =0. (64) 

Our result is that {Ui + U2){I + U\_) should have no anti-Hermitean traceless part. 
Note that this condition has the correct formal continuum limit: writing Ui^2 as 
exp(za_LQi^2) and as exp(ia^a^), we have, for small a_\_, the result a_\_ = qi + 02, 
as required. The above condition in Eq. (|6^ ) can easily be resolved in the SU(2) 
case but we have not yet found a simple closed form expression for Nc > 2. For 
SU(2), one obtains for the initial condition 

U^ = {Ui + U2){ul + ul)-\ (65) 

For the A- field, the singularities arise from the Abelian part of the term 
in the action whose variation with respect to A'^'"' gives 

^Tra^d+{d-A+ - d+A^) . (66) 

Its most singular part is aj9{x^)6{x^). Varying the ±,_L terms (Eq. ^9|) in the 
action (Eq. ^8|) with respect to A'^'"' and selecting the contributions containing 
derivatives, one obtains eventually the result 

"7 = 7^ E Tr^7 {[(Ui - U2){W -I)- h.c.],,„ - [(t/t - I)iUi - U2) - h.c.],-„„, 
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It is easily seen that the above equation has the correct formal continuum limit. 
Writing again Ui^2 as exp(m_Lai^2) and U as ex.p{ia±a±), one finds in the limit of 
smooth fields, a = iJ2n['^iT'^2]n, as required. 

The lattice action is essential to obtain the initial conditions for the evolution of 
fields in the forward light cone. For the evolution, we need the lattice Hamiltonian. 



It is obtained by performing a Legendre transform of Eq. 58 following the standard 
Kogut-Susskind procedure [^]. The analog of the Kogut-Susskind Hamiltonian 
here is 

+ i7 E Tr {^j - U,,n^,^nUl^j ' + J E Tr P|, (68) 

hfi j 

where Ei are generators of right covariant derivatives on the group and Uj^n is a 
component of the usual SU(2) matrices corresponding to a link from the site j in the 
direction h. The first two terms correspond to the contributions to the Hamiltonian 
from the chromoelectric and chromomagnetic field strengths respectively. In the 
last equation <I> = <I)"cj" is the adjoint scalar field with its conjugate momentum 
p = p^-a"". 

Lattice equations of motion follow directly from of Eq. For any dy- 
namical variable v with no explicit time dependence v = {Hl,v}, where v is the 
derivative with respect to r, and {} denote Poisson brackets. We take Ei, Ui, pj, and 
$j as independent dynamical variables, whose only nonvanishing Poisson brackets 
are 

K, = 6,,6ab; {E^, C/^} = -i^imUia''; {Ef, E^} = 26imeabcEf 

(no summing of repeated indices). The equations of motion are consistent with a 
set of local constraints (Gauss' laws). 

The results of this section can be summarized as follows. The four independent 
dynamical variables are Ei, U±, pj and ^j. Their evolution in r after the nuclear 
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collision is determined by Hamilton's equations above and their values at the initial 
time T = are specified by the following initial conditions: 

C/±|r=0 = {Ul + U2){Ul + Ulr^ ■ El\r=0 = • 

Pj\r=o = 2a ; = 0, (69) 



where and a are given by Eq. 65 and Eq. respectively. 



6 Results for gluon production in high energy nuclear 
collisions 

In this section we will discuss recent results for the energy density e as a function 
of the proper time r Work on computing number distributions is in progress 

and will be reported at a later date |^3|. In an earlier work, we confirmed that, 
in weak coupling, the results from our numerical simulations agreed with lattice 
perturbation theory |20|. 



The computation of energy densities on the lattice is straightforward. Our main 



result is contained in Eq. |70|. To obtain this result, we compute the Hamiltonian 
density on the lattice for each p^, and then take the Gaussian average (with the 
weight p?) over between 40 p trajectories for the larger lattices and 160 p trajectories 
for the smallest ones. 

In our numerical simulations, all the relevant physical information is contained 
in g^p and L, and in their dimensionless product g^pL [^^. The strong coupling 
constant g depends on the hard scale of interest; from Eq. ^, we see that p depends 
on the nuclear size, the center of mass energy, and the hard scale of interest; is 
the transverse area of the nucleus. Assuming g = 2 (or as = I/tt), p = 0.5 GeV 
(1.0 GeV) for RHIC (LHC), and L = 11.6 fm for ^n~nuclei, we find g'^pL ^ 120 
for RHIC and ~ 240 for LHC. (The latter number would be smaller for a smaller 
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value of g at the typical LHC momentum scale.) As will be discussed later, these 
values of g^^L correspond to a region in which one expects large non-perturbative 
contributions from a sum to all orders in Qs ~ 6asH, even if as <^ 1. (Recall the 
definition of the saturation scale in Eq. [l^.) Deviations from lattice perturbation 
theory, as a function of increasing g^^L, were observed in our earlier work pO[| . 

We shall now discuss some of the results from our numerical simulations. In 
Fig. 5, we plot ET / {g'^ fi)"^ , as a function of g^fJ-r, in dimensionless units, for the 
smallest, largest, and an intermediate value in the range of g^^Us studied. The 
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Figure 5: er /{g"^ ^)'^ as a function oig^fir for g^fiL = 5.66 (diamonds), 35.36 (pluses) 
and 296.98 (squares). Both axes are in dimensionless units. Note that er = at 
r = for all g^^L. The lines are exponential fits a + /?e~'>"^ including all points 
beyond the peak. 

quantity er has the physical interpretation of the energy density of produced gluons 
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dE jl? jdr] only at late times-when t ^ t. Though et goes to a constant in all three 
cases, the approach to the asymptotic value is different. For the smallest g^^L, et 
increases continuously before saturating at late times. For larger values of g^fJ-L, et 
increases rapidly, develops a transient peak at r ~ and decays exponentially 

there onwards, satisfying the relation a + /?e~'>"^, to a constant value a (equal to 
the lattice dE / L'^ / drjl) . The lines shown in the figure are from an exponential fit 
including all the points past the peak. This behavior is satisfied for all g^fiL > 8.84, 
independently of N. 

Given the excellent exponential fit, one can interpret the decay time td = 
l/^/g'^fj, as the appropriate scale controlling the formation of gluons with a phys- 
ically well defined energy. In other words, td is the "formation time" in the sense 
used by Bjorken In Table 1, we tabulate 7 versus g^^J-L for the largest N x N 
lattices for all but the smallest g^fiL. For large g'^^L, the formation time decreases 
with increasing g^fJ-L, as we expect it should. The reason the smallest value of g^fJ-L 
does not have a transient peak is likely because in this case the kt modes do not 
sufficiently sample the region kf < Qs where non-linearities are important. The few 
modes there are, lie in the perturbative region where the fields can be linearized at 
r = 0. 

The physical energy per unit area per unit rapidity of produced gluons can be 
defined in terms of a function /{g'^fiL) as 

^^ = ^f{g'^L){g'f,f. (70) 
L'^ dr] g'^ 

As discussed in Ref. [21], the function / is obtained for each fixed value of g^fJ.L, by 



taking the continuum limit, i.e., extrapolating g^fia — > 0. In Fig. 6, we plot the 
striking behavior of / with g^fiL. For very small g^fiUs, it changes very slightly but 
then changes rapidly by a factor of two from 0.427 to 0.208 when g^fiL is changed 
from 8.84 to 35.36. From 35.36 to 296.98, nearly an order of magnitude in g^fiL, it 
changes by ~ 25%. The precise values of / and the errors are tabulated in Table 1. 
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Figure 6: eT/{g'^fi)'^ extrapolated to the continuum limit: / as a function of g^fJ-L. 
The error bars are smaller than the plotting symbols. 



What is responsible for the dramatic change in the behavior of / as a function of 
g^fiLl In = gauge, the dynamical evolution of the gauge fields depends entirely 
on the initial conditions, namely, the parton distributions in the wavefunctions of 
the incoming nuclei In the nuclear wavefunction, at small x, non-perturbative, 
albeit weak coupling, effects become important for transverse momenta Qs ~ 6as/i. 
Now on the lattice, pt is defined to be 27rn/L, where n labels the momentum mode. 
The condition that momenta in the wavefunctions of the incoming nuclei have sat- 
urated, Pt ~ GasfJ-, translates roughly into the requirement that g^fiL > 13 for 
n = 1. Thus for g^^iL = 13, one is only beginning to sample those modes. Indeed, 
this is the region in g^iiL in which one sees the rapid change in /. The rapid de- 
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5.66 


8.84 


17.68 


35.36 


70.7 


f 


.436 ± .007 


.427 ± .004 


.323 ± .004 


.208 ± .004 


.200 ± .005 


7 




.101 ± .024 


.232 ± .046 


.165 ± .013 


.275 ± .011 




106.06 


148.49 


212.13 


296.98 




f 


.211 ± .001 


.232 ± .001 


.234 ± .002 


.257 ± .005 




7 


.322 ± .012 


.362 ± .023 


.375 ± .038 


.378 ± .053 





Table 1: The function / = dE/L^/dr] and the relaxation rate 7 = l/r/5/g^/x tab- 
ulated as a function of g^jiL. 7 has no entry for the smallest g'^fiL since there 
eT/{g'^fi)^ vs g^/iT differs qualitatively from the other g^fiL values. 

crease in / is likely because the first non-perturbative corrections are large, and 
have a negative sign relative to the leading term. Understanding the later slow rise 
and apparent saturation with g^fiL requires a better understanding of the number 
and energy distributions with pt. This work is in progress and will be reported on 



separately [53| 



Our results are consistent with an estimate by A. H. Mueller [|57| for the number 
of produced gluons per unit area per unit rapidity. He obtains dN/L?' /drj = c {N^ — 
l)Q'^/'iTT'^ as Nc, and argues that the number c is a non-perturbative constant of 
order unity. If most of the gluons have pt '-^ Qs, then dE/L'^/drj = c' {N^ — 
1) Qg/Ai:'^ as Nc which is of the same form as our Eq. ^ In the g^^iL region of 
interest, our function / ~ 0.23-0.26. Using the appropriate relation between Qs 
and we obtain c' = 4.3-4.9. Since one expects a distribution in momenta 
about Qs, it is very likely that c' is at least a factor of 2 greater than c-thereby 
yielding a number of order unity for c as estimated by Mueller. This coefficient can 
be determined more precisely when we compute the non-perturbative number and 
energy distributions. 

We will now estimate the initial energy per unit rapidity of produced gluons at 
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RHIC and LHC energies. We do so by extrapolating from our SU(2) results to SU(3) 
assuming the Nc dependence to be (A'^^ — 1)/Nc as in Mueller's formula. At late 
times, the energy density is e = (g'^fi)'^ f{g^fiL)j{g'^fj,L)/g'^, where the formation 
time is td = l/^{g'^fiL)/g'^fi as discussed earlier. We find that e^^^^ « 66.49 
GeV/fm'^ and e^^^ « 1315.56 GeV/fm^. Multiplying these numbers by the initial 
volumes at the formation time td, we obtain the classical Yang-Mills estimate for 
the initial energies per unit rapidity Et to be E^^^^ k. 2703 GeV and Ei^^^ k. 
24572 GeV respectively. 

We have compared these numbers to results presented recently [^] for the mini- 
jet energy (computed for pt > Psat, where Psat is a saturation scale akin to Qs)- He 
obtains E^p^^^ = 2500 GeV and Ej.^*-" = 12000. The remarkable closeness between 
our results for RHIC is very likely a coincidence. The Finnish groups results include 



K factor estimates range from 1.5-2.5. If we pick a recent value of iiT ~ 2 [59], we 
obtain as our final estimate, E§"^'^ 5406 GeV and E^^'^ 49144 GeV. 

To summarize, we discussed in this section a non-perturbative, numerical com- 
putation, for a SU(2) gauge theory, of the initial energy, per unit rapidity, of gluons 
produced in very high energy nuclear collisions. Extrapolating our results to SU(3), 
we estimated the initial energy per unit rapidity at RHIC and LHC. We plan to 
improve our estimates by performing our numerical analysis for SU(3). Moreover, 
computations in progress to determine the energy and number distributions should 
enable us to match our results at large transverse momenta to mini-jet calcula- 
tions S. 



7 Summary 

In these lectures, we have discussed a classical effective field theory approach to 
scattering at very high energies. At these energies, a saturation scale Qs{x) controls 
the dynamics of high energy scattering. How this scale changes as we go to small x 
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is described by a non-linear renormalization group equation |15, The solutions 
of the RG equations and the inclusion of effects such as the running of the coupling 
in the regime of strong non-linear fields need to be better understood. In particular, 
one would like to investigate possibly striking experimental signatures of this regime. 
Since Qs{x) » Aqcd, weak coupling methods may be applicable. In these lectures 
we have shown how one may begin to apply these weak coupling methods to study 
DIS and high energy scattering. 
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